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Abstract: The classical differential cross-section is calculated for high-energy small-angle 
gravitational scattering in the factorizable model with toroidal extra dimensions. The 
three main features of the classical computation are: (a) It involves summation over the 
infinite Kaluza-Klein towers but, contrary to the Born amplitude, it is finite with no need 
of an ultraviolet cutoff, (b) It is shown to correspond to the non-perturbative saddle- 
point approximation of the eikonal amplitude, obtained by the summation of an infinite 
number of ladder graphs of the quantum theory, (c) In the absence of extra dimensions it 
reproduces all previously known results. 
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1. Introduction 

The search for large extra dimensions (LED), especially in view of the forthcoming experi- 
ments at LHC, constitutes a very exciting direction in high energy physics. After precursor 
ideas about the Universe as a topological defect in higher-dimensional space-time |^ and 
the proposal of TeV-scale internal dimensions related to supersymmetry breaking in string 
theory large extra dimensions were discussed in several contexts. Among the concep- 
tually and technically simplest is the ADD picture |^], according to which the standard 
model particles live in a four-dimensional space (the brane), embedded in a L)-dimensional 
bulk inhabited only by gravity, and with the extra 6 = — 4 dimensions compactified on 
a torus. The D-dimensional Planck mass M^, is supposed to lie in the TeV region and the 
LED have submillimeter size. In this scenario there is an infinite tower of massive Kaluza- 
Klein (KK) gravitons whose existence may be detected at present and future colliders 
Q , either as missing energy in collisions due to emission of KK gravitons (which are weakly 
interacting after being created), or via processes which would otherwise be impossible or 
very much suppressed. 

Quantum collision processes with exchange of virtual KK gravitons are a very useful 
tool to test the model [^. Unfortunately, tree diagrams containing the propagator of the 
KK tower diverge in the ultraviolet (UV) due to an infinite sum over the KK graviton 
masses The divergence appears already at tree level and is due to the emission of 
infinite momenta into the compactified dimensions. Several proposals were made in order 
to cure this problem. One was to introduce a UV cutoff |Q, ^, another to take into 
account brane oscillations and the associated Nambu-Goldstone modes [^, another yet to 
introduce brane thickness |§, or finally, to sum up ladder diagrams within the eikonal 
approximation Q. However, none of them solves the problem completely. The UV cutoff 
at Planckian energies while it certainly does exist, does not look very natural in this 
context. It would mean that we are able to probe the Planckian regime using (relatively) 
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low energy processes. Appealing to features associated with the physical brane, such as 
tension and thickness, introduces additional elements into the ADD model, which seem 
to spoil its overall consistency. Indeed, for a brane of non-zero tension one can not use 
any more the flat Minkowski background. Instead, one has to deal with solutions of the 
Einstein equations, which are non-trivial even in the absence of matter on the brane. In 
the case of codimension one, such a solution is well-known and is the basis of the Randall- 
Sundrum II model [^], with a different graviton spectrum. Finally, the eikonal calculation 
amounts to using the Fourier-transform of the Born amplitude and thus, it also suffers 
from ambiguities associated with the divergence of the latter. Indeed, the evaluation of the 
eikonal phase gives a finite result if a certain order of integration is used in calculating the 
Fourier integral, while it is divergent if one merely changes the order of integration. 

On the other hand, it has been argued that at ultrahigh energies, particle scattering 
in four dimensions not only becomes dominated by gravity, but in addition it involves only 



classical gravitational dynamics 12, Indeed, quantum gravity effects should 

not, by definition, be important in the classical limit h ^ 0. This, in terms of the two 
relevant lengths, i.e. the Planck length /pi = {hGi/c^Y/"^ = h/Mpic and the gravitational 
radius associated with the energy of the collision = GA^^fsjc^^ implies that Vg ^> /pi, 
which is equivalent to the condition ^ My-\(? of transplanckian energies. Thus, to 
study the scattering of two point-like particles with Planck energy, 't Hooft used a shock 
wave approximation for the field of the moving particle and obtained a result similar to the 
Veneziano amplitude. Later on, it was shown that in four dimensional quantum gravity 



1 15], as well as in string theory ||T^, the eikonal approximation (free in both cases of the 
ambiguities mentioned above) reproduced the result of 't Hooft. 

The Planck length and the gravitational radius r* in the D-dimensional ADD model 
are, correspondingly 

where is the TeV-scale mass parameter. The above reasoning remains essentially the 
same and shows that in the transplanckian regime ^/s ^ M^^c^ scattering is also classical, at 
least for some range of momentum transfer. Moreover, the particularly exciting proposals 
of black hole creation at LHC or in cosmic rays [|^] were based on a purely classical 
picture. 

Although the eikonal approximation of particle scattering in ADD has been discussed 
by a number of authors |^], |^], no classical calculation of the cross-section was found in 
the literature. The purpose of the present paper is to fill in this gap, and in addition to 
provide an independent check of the validity of the eikonal approach, whose applicabil- 
ity in the context of ADD is not yet rigorously proven. Indeed, it is hereby demonstrated 
explicitly that the classical theory reproduces the saddle-point result of the eikonal approx- 
imation and is essentially non-perturbative in the quantum sense. Thus, classical theory 
gives non-trivial results in ADD, in contrast to four dimensional gravity where the classical 
elastic scattering cross-section coincides with the Born approximation and is perturbative. 
Finally, it should be pointed out that the reason the classical cross-section is finite, even 
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though it too involves an infinite summation over massive KK modes, is because in the 
classical calculation the potentially divergent integrals contain oscillating factors, which 
effectively cutoff the modes that cannot be excited by the source. This applies to all kine- 
matical regimes and, in particular, is analogous to the situation in the classical derivation 
of Newton's law [|19|. 



2. The setup 

Consider the Fierz-Pauli lagrangian in D-dimensional Minkowski space, with 5 = D — 4 of 
its spatial dimensions being a torus with equal radii R 

c = -h^^^nhMN + \hnh - h'^^dMdNh + h^^^dMdph^ - ^h'^^Tmn. (2.1) 

where M, .. = 0, 1, 2, D — 1. The Minkowski metric is tjmn = diag(l, —1, —1, ...) and 
□ = ■q'^^^ dudN ■ The gravitational field huN is coupled to a conserved matter stress-tensor 
Tmn {OnT^'^ = 0). All fields {h = n^^^huN = h%) are functions of 

x^^ = (x^y*), ^ = 0,...,3, i = l,...<5, (2.2) 

and are periodic under the translations yj — > yj + 27ri?, j = 1, . . . , 5. Thus, for instance 

ni=— oo n^=— oo ^ ^ ^ 

where V = {2'kRY is the volume of the torus. In what follows we abbreviate the sum over 
the KK modes as and denote the momenta transversal to the brane as p^rp = v} jR. 
The tensor hMN{x^) has been split into an infinite sum of four-dimensional KK modes 
h^j^{x^^) with (mass)^ equal to p^. 

-¥ 



In the harmonic gauge d^h^^^ = \d^h the Einstein equations for Hmn read: 



□/lAfAT = — {tmn — ^ _^ = —Jmn, T = r]'^^ Tmn- (2.4) 

According to the ADD scenario it is assumed that the matter stress-tensor is localized on 
the brane, located at y = 0, and carries only four-dimensional indices: 

TMiv(x^) = vtfVNT^.u{x)S'iy). (2.5) 

It is worth noting, however, that this assumption is consistent with Einstein's equations 
only at the linearized level, with matter dynamics governed by the non-gravitational forces 
alone. Given ( |2.5D it is consistent to set the graviphotons hi^ and the non-diagonal part 
of the scalar matrix hij to zero. Finally, the diagonal components of hij are all equal and 
generated by the trace of the energy- momentum tensor: 

T 

Dhij = -^Dj^^ij- (2.6) 
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Their zero modes = are the so called radions, which describe deformations of the 
torus caused by the presence of matter on the brane. The massive components 7^ of 
^MAT could be rearranged into the massive four-dimensional graviton and massive scalars 
in the usual Higgs mechanism language |Q, Q, but this is not necessary for the present 
discussion. 

The Z)-dimensional Planck mass M^, is defined by >f|j = IGvr/M^^^ = IGttGd and is 
related to the four-dimensional one Mpi via 

The retarded Green's function of the D'Alembert equation satisfies 

OGDix, x', y - y') = -S\x - x')6^{y - y'). (2.7) 

Its Fourier transform reads: 

Oo(. - y - y-) = / Ae-^ E /"^^l^ ■ P.B) 

The solution of ( p.^ ) with the source localized on the brane 

JMN{x,y) = JMN{x)S^{y), (2.9) 



IS 



hMN{x,y) = j Gd{x - x' ,y - y')JMN{x' ,y')d'^x'dy' 



{2ttYV 



E 



g«PT-y 

2 _ 2 I ■ "^MNix') d x'. 



Its restriction to the brane 



hMN{x) = hMN{x,y) 
can be rewritten using the amputated propagator 



y=0 



(2.10) 



(2.11) 



hMN{x) = J G4(x - x')JMN{x')d'^x', G^ix - x') = Gd{x -x',y - y') 
The momentum space four-dimensional retarded propagator thus reads 



y=y'=o 



• (2.12) 



G^ip) = 77 ^ : — f7- 

^ P Pt ~^ ^^P 



(2.13) 



Equivalently, taking the four-dimensional Fourier transform, defined by "^{k) = J '^{x)e^^'^d'^x 
and \I'(x) = f ^ {k)e~^^'^ d'^k / {2tt)^ of the above equations, the retarded solution of Eq. (^ 
becomes 

hMN{k) = —2_^ ji^ 2 , , ,n • (2-14) 



V 



k"^ — p^ + iek^ 
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3. The ultra-relativistic elastic scattering cross-section 

Consider next the small angle ultrarelativistic scattering of two particles on the brane, with 
masses ra and m' respectively, interacting via D-dimensional gravity. Using for notational 
simplicity the same parameter r for both trajectories, a convenient way to describe them 
is 

z^^[T) = z^ + ?^T + 5z^'{T), ^'A*(r) = + + fc'^(r), (3.1) 
m m' 

with 5z, 5z' treated perturbatively. The asymptotic values of their momenta are 

P^ = m lim i^(r)=p^+m lim 5i^(r), P'^' = m' lim i'^(r) = p'^ + m' lim di'^'ir), 

P'' = m lim i^(T) = p^' + m lim 5i^(r), P'^' = m lim i'^(T) = p^' + m' lim 5z^{t), 

(3.2) 

while momentum conservation implies P'^ + P'^ = P^ + P'^. The Mandelstam variables 
s, t are 

2 



s = {P + P')\ t = q' = {P-Py = m^ [ lim di^'ir)- lim (5i^(r) , (3.3) 

\ T^oo T^— oo y 

and we consider the ultrarelativistic regime and the small-angle approximation 

s > m^ m'^ (3.4) 

in which case s = {p + p')"^. 

Since the momentum transfer depends only on the deviation 5z^ of one of the 
particles, it will be convenient to work in the rest frame of m' before collision. In that 
frame 

p'^ = m'(l, 0,0,0), p^ = m7(l,0,0,^;), 7 = —!=, (3.5) 

V 1 — "V 

and with no loss of generality one may set in addition 

4^ = 0, ^0^ = = (0,6,0,0), (3.6) 

where h is the impact parameter. 

The particles equations of motion following from the (^]5) are D-dimensional geodesic 
equations in the metric 

Qmn = r]MN + >(DhMN, (3.7) 

but it is easy to show that the particle moving on the brane in zero order in >co will remain 
on the brane. So the matter stress-tensor of the two particle system is four-dimensional: 



j [mz^'z''5^{x - z{t)) + m'i'^'i"'5\x - z'{t))] dr. (3.8) 



Leaving aside the classical mass renormalization needed to take into account the grav- 
itational self-action, we have to take as h^u in the equation for the particle m the retarded 
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field of the partner particle m' . In the small-angle approximation it is assumed that the 
deviation from the unperturbed rectilinear motion is small, and we get perturbatively 

n^-ji;, = x,,n^- (^/i,A,p - ^/iAp,.) (3.9) 

where Yi^'^ is the projector onto the space transverse to the momentum 

= r//^^ - (3.10) 

The gravitational field h^^ of the particle m' is given by ( p.l4| ), where we have to 
substitute the Fourier transform of the second term of the stress-tensor ( |3.8D 

Using the reparametrization invariance of the particle trajectories one may choose r in 
such a way that p^5z^ = for both particles. In this gauge the solution of (|3.9|) becomes 

= ^ / ,'!^'^f7"^:ij rA. (3.12) 

The vector has the form 

Q/^ = Ak^ + Bp^" + Cpl^, (3.13) 

where 

Upon differentiation with respect to r and integration over using that in the chosen 
Lorentz frame b(k ■ p') = 6{m'k^), one obtains 



9 f ik-b—ik-pr /m 

Sz^'ir) = / y 



Sk. (3.15) 

fcO=0 

To calculate the momentum transfer in ( p. 3D we need the asymptotic values of 5z^^ as 
r iboo, which one computes next. In the chosen Lorentz frame the exponent in ( ^.15 ) is 

k ■ b — k ■ pr/m = kxh — k^vyr. (3.16) 

(a) Start with the integral over kz in ( p. 15 ). Define 

^,2^ 2 dh, >.^ = kl + kl+pl. (3.17) 

-oo + X 

The terms B and C, as well as the one proportional to Ak^ vanish in the limit r ±00, 
because 

/ ^dkz = -e-^'^l^l ^ , as r ^ ±00. (3.18) 

7_oo ki + >i:^ K 
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The term proportional to Akg is zero, because it is evaluated at ko = 0. The y component 
vanishes by parity. So, the only component left is the one proportional to Akx- In this 
case, parity implies that only the sine part of the exponential contributes. Use 

hm ""i^-^""^ = ±n5ik.jv), (3.19) 



to perform the kz integration. Then, insert A from (3.14) to obtain 

hm F(t) = ±-^^ ( {p ■ p'f - ■ (3-20) 



(b) Next, denote by = + Pr and perform the integral over k 



dkx = 7^e"■^^ (3.21) 



oo ~l~ 



(c) Insert this into (3.15) and replace the sum over KK modes by a continuous integration 

VSs^i r 5-1, n 27r^/' 



(6/2) ' 



(3.22) 



to obtain 

di"(±oo) = ± — ^ ( (p ■ p')' - r dky r e-^^ p^r'dpT- 

^ ' 22+^7r^/2m'm27i;r(5/2)M2+^ V 5 + 2 ) J _^ " Jq 

(3.23) 

(d) To perform the remaining integrations pass to polar coordinates pT = i^cosa, ky = 
X sina, p^jT^dkydpT = K^dK cos^~^ ada, and integrate over K from zero to infinity and 
over a from —tt/2 to tt/2. The asymptotic values of the velocity of m in the field of m' are 
then: 

x-x,^ ^ ^ 2r(5/2 + l) f, ,.2 m^m'^\ ,^ 
5z'' (+00) = +— — --^ — —J (p ■ p'Y ; . (3.24) 

Inserting this into ( p. 3D we find for the square of the momentum transfer 



-t - 2''T^(5/2 + 11 ^ f ^P-P'^^ - J—V (3 25) 

In the ultrarelativistic limit 7 ^ 1 , v ~ 1 , one has s = 2p ■ p' ^ mm' and the above 
expression simplifies to: 



_ 22r2(,5/2 + 1) 
~* ~ 7r^(M,6)2{5+i) M| ■ 

Finally, the differential cross-section, defined as usual by da = 27rbdb, is 



(3.26) 



da 1 / s2 X 1/(5+1) 



In particular, for 5 = Q 

dt 



da _ 47rG|s2 



(3.28) 
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which coincides with the well-known formula for small angle scattering of General Relativity 

The scattering angle 9 is given by tan^ = -v/— t/m7u, thus, small scattering angles 
mean \t\ <C m'^j'^v'^, which for ultrarelativistic velocities gives the range of validity of our 
approximation 

{M.bY+' » ^. (3.29) 
4. Relation to the eikonal approximation 

Consider the elastic scattering of massive scalar particles on the brane in the high-energy 
limit s ^ m^. The Born amplitude contains the t-channel propagator involving the sum 
over Kaluza-Klein modes. Passing to the continuous integration over the momentum of 
gravitons pt in extra dimensions we have 0]: 



2(2^7 -t+py 



MBorn{s,t) = / " , . (4.1) 



This integral in the general case requires a UV cutoff. An alternative way to get the final 
amplitude for the small-angle high-energy scattering is to use the eikonalized form of the 
amplitude |, |, g 

7Weik(s, t) = 2is j e''!-^ (l - e^'^^^'^)) d%, (4.2) 

where the two-dimensional vectors q, b lie in the transverse plane, with b the impact 
parameter vector. The transverse component q of the momentum transfer in this approxi- 
mation satisfies q^ ~ —Q'^Qin so that t ~ — q^- This expression in the usual four-dimensional 
theory corresponds to summation of the ladder and crossed-ladder diagrams (for a detailed 
calculation within the quantized linearized General Relativity see fl^). Actually, this in- 
volves UV divergent loop diagrams, but it can be shown that the leading contribution 
in the high-energy limit is independent on the cutoff. In the ADD linearized gravity the 
situation is believed to be the same, though no explicit analysis is available. Therefore 
our classical calculation provides an independent check of the applicability of the eikonal 
approximation in the ADD framework. 

The Born amplitude corresponds to the first term in the expansion of the exponential 
mil) 

MBornis, t) = 2s j e'^-\{s, b)d\ (4.3) 
and is used to extract the eikonal x its inverse Fourier-transform 

X{s, ^) = Ysj ^^''''''^^oUs, t) (4.4) 



Notice, that although the Born amplitude itself may be divergent, the integral (4.4) is finite 
if one first integrates over q, but not px- Indeed, choose the coordinates as in the previous 
section to write 
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and integrate first over qx using a contour integration which gives an exponential factor 
cutting the potentially divergent integral over px- Then, integrate over qy = xcosa and 
Pt = x: sin a to obtain 



X{s,b) = / dxj dae-'^'^Ss-icos^-^ax^ = i^jj , (4.6) 



where 



Then, the eikonal amplitude ( |4.2| ) becomes 

Meikis, t) = 4ms j Joiqb) (l - e^'^^''^)) bdb. (4.8) 

The unity in the parenthesis gives no contribution. In the remaining part and in the regime 
qbc S> 1 of interest here, one may replace the Bessel function by its asymptotic to obtain 

/•OO Oh r 1 
Meik{s,t) =27ris JL ^iMQ,b)-i7r/A^^ii^-{q,b)+tn/4 (4^9) 

Jo V 

where ^ 

^l^^(q,b)=±qb+(^) . (4.10) 



This can be evaluated using the stationary phase method. The exponent of the second 
term having no critical points of first order in the domain of integration, is ignored. The 
stationary point for the first exponent is 

^ - b-('Ji] r ib ) - (4 11) 

db b,~^' V <? y ' ^ ~ 6f+' ' 

and upon integration 



s 



7Weik(s,t) = , e^^^- = -^e^^''-. (4.12) 

^ ' ,J6{6 + 1) bi'^ V5Tl ^ ^ ' 



jS+l 



Substituting be from (|4.7|) one obtains 

4V^se^('?^'--'^/2) ^x2,sr((5/2 + 1)Y^ _ 40Fse*(''''=-'^/2) ^2^sr(5/2 + 1)Y^ 



qV6TT V 80Fg J qV6TT V M^+^q 

(4.13) 

Note, that although a few intermediate steps are singular for 5 = and seem to require a 
separate discussion, the final formula is valid for (5 = as well. 
The corresponding cross-section reads 

As advertised, it is identical to our classical result. 
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5. Conclusions 



A purely classical calculation was presented of the high energy elastic scattering cross 
section in the ADD scenario. Our approach is entirely free of the ambiguities associated 
with the divergence of the Born amplitude with the virtual graviton exchange typical for 
ADD. Ultrarelativistic small angle gravitational collision in four dimensions is a special 
case, in which it agrees with 't Hooft' s result, which in turn coincides with the Born 
quantum cross-section. In the presence of extra dimensions it was shown that the lowest 
order small angle classical approximation reproduces the essentially non-perturbative result 
of the quantum eikonal calculation in the saddle-point approximation. Thus, the classical 
computation in the above kinematical regime is non-trivial, unambiguous as well as reliable 
and, therefore, worth applying to other processes like bremsstrahlung |22]. 
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